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Abstract 



Data collected at centre-of-mass energies of 181-184 GeV by ALEPH at LEP, 
. corresponding to an integrated luminosity of 56.9 pb~^, are analysed in a search 

r^ I for pair-produced charged Higgs bosons H^. Three analyses are employed to select 

C^ ' the T^i'rT^i^T, csT^Vr/csT^Vr and esse final states. No evidence for a signal is found. 

Mass limits are set as a function of the branching fraction B(H^ — ^r^i/,-). Under the 
assumption that the decay modes considered cover the totality of the possible final 
states, charged Higgs bosons with masses below 59 GeV/c^ are excluded at 95% C.L. 
independently of B(H+ -^t^Vt). 
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1 Introduction 

Despite the success of the Standard Model of electroweak interactions in describing 
experimental observations, not much information is available about its cornerstone, the 
Higgs sector. In its minimal version, the Higgs mechanism is implemented by adding 
only one doublet of complex scalar fields, resulting in one additional physical scalar 
state, electrically neutral, commonly referred to as the standard Higgs boson. The most 
important phenomenological consequence of an extended Higgs structure is the appearance 
of additional physical spin-0 states ^ . For example, with the addition of one more doublet 
of complex scalar fields, five physical states remain after spontaneous breaking of the 
SU(2)l X U(1)y symmetry to give mass to W^ and Z gauge bosons: three neutral and a 
pair of charged bosons. Among the possible choices, multi-doublet models are theoretically 
interesting because they automatically lead, at tree level, to m-w = rnz cos 6'w and to the 
absence of flavour changing neutral currents, two major constraints which must be satisfied 
by any extension of the Standard Model to agree with the experimental observations. 

This letter describes a search for pair production in e^e~ collisions of the charged Higgs 
bosons H^ predicted in two-Higgs-doublet extensions of the Standard Model. The analysis 
uses the total integrated luminosity of 56.9 pb^^ collected in 1997 with the ALEPH detector 
at LEP, at centre-of-mass energies from 181 to 184 GeV, hereafter called the 183 GeV data. 

Pair production of charged Higgs bosons occurs mainly via s-channel exchange of a 
photon or a Z boson; in two-doublet models, the couplings are completely specified in 
terms of the electric charge and 6^^, making the production cross section depend only 
on one additional parameter, the charged Higgs boson mass ?71h±- As expected in most 
implementations of multi-doublet models |]I[ , it is assumed that H"*" decays, with negligible 
lifetime, predominantly into cs or t^Ut- (and the respective charge conjugates for H^). 
Additional decay channels, such as those involving neutral Higgs bosons, are not considered 
here. Since the relative weight of the two main channels depends on the details of the model, 
no assumption is made about the decay branching fractions, and three different selections 
are developed to address the possible final states esse, csr'Ur/csT'^i'r (hereafter referred to 
as csr~z/T-) and t+z/t-t"//,-. Under the same hypothesis, the negative results of the searches 
performed using 27.5 pb~^ collected at centre-of-mass energies ranging from 130 to 172 GeV 
allowed ALEPH to exclude charged Higgs boson masses less than 52 GeV/c^ at 95% C.L., 
independently of the final state 0. Using the data recorded at the same centre-of-mass 
energies, an excluded domain up to 54.5 GeV/c^ has also been reported by DELPHI 0. 
Charged Higgs boson masses up to 57.5 GeV/c^ and 59.5 GeV/c^ have been excluded by 
L3 [^ and OPAL P , respectively, using their data recorded at centre-of-mass energies up to 
183 GeV. Less general limits have also been set by ALEPH @], CLEO and CDF [g. 

The letter is organized as follows. After the description of the relevant parts of the 
ALEPH detector in Section 2, the event selections are detailed in Section 3. The results and 
the conclusions are given in Sections 4 and 5. 



2 The ALEPH Detector 

The ALEPH detector is described in detail in Ref. [H. An account of the performance of the 



detector and a description of the standard analysis algorithms can be found in Ref. |T0 
Here, only a brief description of the detector components and of the algorithms relevant for 
this analysis is given. 

In ALEPH, the trajectories of charged particles are measured with a silicon vertex 
detector, a cylindrical drift chamber, and a large time projection chamber. These are 
immersed in a 1.5 T axial field provided by a superconducting solenoidal coil. The 
electromagnetic calorimeter, placed between the tracking system and the coil, is a highly 
segmented sampling calorimeter which is used to identify electrons and photons and to 
measure their energies. The luminosity monitors extend the calorimetric coverage down to 
34 mrad from the beam axis. The hadron calorimeter consists of the iron return yoke of the 
magnet instrumented with streamer tubes. It provides a measurement of hadronic energy 
and, together with the external muon chambers, muon identification. 

The calorimetry and tracking information are combined in an energy flow algorithm, 
classifying a set of energy flow "particles" as photons, neutral hadrons and charged particles. 
Hereafter, charged particle tracks reconstructed with at least four hits in the TPC, and 
originating from within a cylinder of length 20 cm and radius 2 cm coaxial with the beam 
and centred at the nominal collision point, are referred to as good tracks. 



3 Event selections 

In order to ensure a good discovery potential independent of the branching fraction 
B(H"'"-^r+z/^), three selection procedures are designed for the topologies r~''z/^r~z/^, csr^z/,- 
and esse. As in Ref. , the most relevant selection criteria for the three selections are chosen 
in order to achieve, on average and in case no signal is present, the best 95% CL. limit 
on the H"'"H~ production cross section. To do so, each selection is optimized individually 
with the most optimistic B(H"'"— >-r"'"^'T-) in each case (100%, 50% and 0% for the T^Vr'T~i'T) 
cst~Ut- and esse channels, respectively, for which the combined contribution of the other 



two analyses is minimal), following the prescription of Ref. ITTI modified to include the 
possibility of partial or full background subtraction. In the following sections, the Monte 
Carlo samples used in designing the selections are described and the changes with respect to 
the analyses published in Ref. are presented. In each case, the subtractible background, 
i.e., that for which the theoretical knowledge and the simulation accuracy are considered 
to be under control, is estimated together with the related systematic uncertainty. 



3.1 Monte Carlo Samples 

Fully simulated Monte Carlo event samples reconstructed with the same program as the 
data have been used for background estimates, design of selections and cut optimization. 



Samples of all background sources corresponding to at least 20 times the collected luminosity 
were generated. The most important background sources are e"'"e^ -^ t^t~ , qq and four- 
fermion processes (including W^W~production), simulated with KORALZ [12], PYTHIA [13] 



and KORALW |T4] 



The signal Monte Carlo events were generated using the HZHA ][T^ generator, extended 



for charged Higgs boson production as described in Ref. ^. Samples of at least 1000 signal 
events were simulated for each of the various final states for charged Higgs boson masses 
between 40 and 80 GeV/c^. 

3.2 The T^i/^T~i>T. final state 

The final state produced by leptonic decays of both charged Higgs bosons consists of 
two acoplanar r's and missing energy carried away by the neutrinos. Since this topology is 
very similar to that of stau pair production, the selection described in Ref. ]]T6[ is used here 



to search for charged Higgs bosons in the T~^UrT~i'r channel. This selection exploits the fact 
that the signal events contain at least four neutrinos, leading to large missing energy and a 
large acoplanarity of the visible system. Background from W"'"W~ production followed by 
leptonic W decays is suppressed by vetoing events with energetic electrons or muons, which 
are softer when originating from r decays. 

Efficiencies to select events from H'^'H" -^ t^v^t^v^ are of the order of 45%, as shown 
in Table || for a representative set of Higgs boson masses. The total background expected 
amounts to 6.5 events, consisting mainly of irreducible background from W^W^ -^ 
T^UrT^Ur- In the data, four events were selected, in good agreement with the Standard 
Model expectation. For the interpretation of the negative result of the search in terms 
of mass limits, the part of the expected background coming from W"'"W~ is subtracted; 
the latter amounts to 5.0 events, including a reduction of 3% to account for systematic 
uncertainties ]|16 . 



3.3 The csT i^r final state 

The mixed final state, csr"//,-, is characterized by two jets originating from the hadronic 
decay of one of the charged Higgs bosons and a thin r jet plus missing energy due to the 
neutrinos from the decay of the other. 

Two complementary approaches are used to select the mixed final state: in one 
selection, called the global analysis, global quantities such as acoplanarity, thrust, and 
missing momentum are predominantly used whereas the second selection, referred to as the 
topological analysis, relies more on the specific r jet reconstruction. As the analyses are 
described in detail in Ref. ]@] , the focus here is on changes other than a simple rescaling of 
the cuts with y/s. 

As wih± approaches ?nw the sensitivity of both analyses is limited by the W"'"W~ 
background. Two changes are introduced to improve the rejection of this background. 
The first change concerns the momentum of the leading lepton (electron or muon) which is 



Table 1: Efficiencies e (in %), numbers of Standard Model background events expected (-/V^^g) 
and subtracted (-/V^w), and numbers of observed candidates (A'obs) for the three analyses at the 
centre-of-mass energy of 183 GeV, as functions of the charged Higgs boson masses (in GeV/c^). 
For the mixed and four-jet channels, numbers are quoted within the windows defined by the sliding 
cuts, therefore implying some overlap among different charged Higgs boson mass hypotheses. 
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now required to be less than 24 GeV/c in both analyses. 

The second change involves the rejection of the ri/qq' final state of W pair events: cuts 
are introduced which depend on the signal mass hypothesis [sliding cuts) and are tightened 
with increasing mH±. 



In the global analysis the rejection is achieved with cuts on the acoUinearity angle r|^f^ 
of the two hadronic jets and their invariant mass m^^. The acoUinearity angle is required 
to be less than [35+mH±/(GeV/c^)] degrees. The event is rejected if the invariant mass lies 
outside the range of [mH± — 10 GeV/c^,mH±+5 GeV/c^]. The good agreement of data and 
Monte Carlo in these two variables is shown in Fig. |I[ 

In the topological analysis the invariant mass mlf^ of the two jets assigned to the 
Higgs boson must lie in a window between mH± — 15 GeV/c^ and r7iH±+5 GeV/c^. Their 
acoUinearity r]H must be less than [40+mH±/(GeV/c^)] degrees. 

The complete set of cuts is listed in Table ||. As in Ref. 0, events are accepted if 
they pass either analysis. Typical efficiencies and background expectations are given in 
Table |l|. The main contributions to the systematic error (3%) on the efficiency are: Monte 
Carlo statistics; the luminosity measurement accuracy (<1%); the uncertainty (<2%) on 
the knowledge of the inefficiency introduced by beam related energy deposits at polar angles 
below 12°, studied using events triggered at random beam crossings. 

Nine events were selected in the data for Higgs boson masses from 40 to 75 GeV/c^, in 
agreement with the background expectation of 8.7. 

For the determination of the result, the W"^W~ background, representing 93% of the 
total contamination, is reduced by a systematic error of 13% and subtracted. Here the 
systematic error represents the statistical precision of a test of the W"'"W~ Monte Carlo 
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Figure 1: Distribution of (a) the invariant mass of the quark jets and (b) their acolUnearity 
angle used in the csT~v'r global selection. The dots are the 183 GeV data, while the shaded 
histogram is the background expectation, normalized to the recorded luminosity. The dashed line 
is the distribution for a signal with a Higgs boson mass of 55 GeV/c^, arbitrarily normalized. 
The arrows indicate the position of the cuts applied for this specific choice of the charged Higgs 
boson mass hypothesis (see text for details). Only a subset of the cuts is applied here to preserve 
sufficient statistics. 



Table 2: Summary of the cuts apphed in the csr Ur analyses. Units for masses and momenta 
are GeV/c^ and GeV/c, respectively. Variables not defined in the text have the same meaning as 
in Ref. [2]. 
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for masses reconstructed in the region of the expected sensitivity. The test is performed 
using a data sample dominated by the W"'"W~ process, selected by relaxing the cuts on the 
hadronic acoUinearity and on the energy of the leading lepton and requiring the hadronic 
mass to be less than 75 GeV/c^. In the data, 63 events were observed, in agreement with 
the expectation of 73.5 events. 

3.4 The esse final state 

For this channel, the hadronic decays of the two charged Higgs bosons lead to a final 
state with four well separated jets. With respect to Ref. P|, the preselection applied to 
identify four-jet final states and the choice of jet pairing are unchanged, but the variables 
discriminating signal and Standard Model processes are exploited in a different way to face 
the larger W"''W~ background. These variables are 

• the production polar angle ^prod between the Higgs boson momentum direction and 
the beam axis; 

• the decay angles 6dec,i (^=1)2) in the rest frame of the two reconstructed H^ 
candidates; 



the chi squared xlc ^^ ^^^ 5C-fit. 



In addition, a c-jet tagging variable Ctag is introduced to take advantage of the presence 
of c quarks in the signal decay products. This variable is the output of a neural network 
trained to discriminate c-jets from light quark jets. The lifetime and specific decay modes 
of D mesons as well as jet-shape properties are exploited. A detailed description can be 



found in Ref. [0 



The four variables are combined linearly into one discriminant observable: 
D = - cos^ ^prod + 0.4 Ctag - 0.2 Min(cos^dcc,*)^ - 0.6 xlc- 
The distribution of D is shown in Fig. |^. Events are accepted ii D > —0.4. 

The selection cuts are summarized in Table ^. After these cuts, the total background 
expected for reconstructed charged Higgs boson masses lower than 75 GeV/c^ amounts to 
42.8 events. Efficiencies are of the order of 35% within a dijet mass window of ±3 GeV/c^ 
around the Higgs boson mass hypothesis, as shown in Table [^. 

As in Ref. , given the level of irreducible background, the sensitivity of the analysis 
is considerably increased by subtracting the expected background from Standard Model 
processes. For this purpose the dijet mass distribution as obtained from the background 
Monte Carlo is parametrized by the sum of a polynomial and a Breit-Wigner distributions. 
The comparison with the data (Fig. |^) shows that the parametrization is consistent with 
the observation. In the following, the subtracted background is conservatively reduced by 
11%, corresponding to the statistical uncertainty of this comparison. 

The systematic error on the number of signal events expected is estimated to be 2%, 
dominated by the Monte Carlo statistical uncertainty, with small additional contributions 
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Figure 2: The distribution of the discriminant variable D used in the esse selection. Shown 
are data (points with error bars), Monte Carlo for the expected background sources (cumulative 
histograms, normalized to the recorded luminosity) and the Monte Carlo expectation for a signal 
with mjji = 55 GeV/c^ (dashed histogram, arbitrary normalization). The arrow indicates the 
position of the cut applied. Some cuts have been relaxed to preserve sufficient statistics. 



from the luminosity measurement and possible inaccuracies in the simulation of the energy 
flow reconstruction. 



Table 3: Summary of the cuts applied in the esse analysis. Units for masses and momenta are 
GeV/c^ and GeV/c, respectively. The preselection variables have been defined in Ref. [2]. 



Four-jet preselection 


Good Tracks 


>7 


Charged Energy 


> io%yi 


IP ■ 1 


< 1.5 mvis-90 


TTTem 


< 90%Ejet 


Y34 


> 0.003 


Thrust 


<0.9 


Equal Mass and Spin-0 constraints 


C'prod; "dcc,i, XsC? '^tag 


D > -0.4 
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Figure 3: Distribution of the dijet mass iii||^ obtained after applying all cuts of the esse selection; 
shown are 130-183 GeV data (dots with error bars), background Monte Carlo (shaded histogram) 
and the background parametrization used in deriving the limits (solid line). 



4 Results 



The number of candidate events observed in the data collected at centre-of-mass energies 
from 181 to 184 GeV are given in Table |l| for different charged Higgs boson masses and 
for each of the three analyses presented in Section 3. A total of 60 events is retained 
for mu± < 75 GeV/c^, consistent with the 57.9 events expected from Standard Model 
processes. Since, in addition, the mass distribution in the esse channel does not show any 
significant accumulation outside the W region (Fig. ^), the results of the three selections 
described in this note are combined with those obtained using 130-172 GeV data to set an 
improved 95% C.L. lower limit on the charged Higgs boson mass, following the procedure 
described in |l^ for the combination of the confidence levels and the prescription of [^ 



for the background subtraction. 

The separate results of the three analyses are displayed in Fig. ^, where the contours 
corresponding to expected (dash-dotted curves) and observed (solid curves) confidence levels 
of 5% (equivalent to a 95% C.L. exclusion) are drawn. For B(H+^r+z/^) = 0, 0.5 and 1, 
values maximizing in turn the weight of the three channels esse, cst~Uj. and t^Ut-t^Ut-, 95% 
C.L. lower limits on rriH± are set to 62, 59.5 and 74.5 GeV/c^. 



The result of the combination of the three analyses is displayed in Fig. ^. Charged Higgs 
bosons with masses less than 59 GeV/c^ are excluded at 95% confidence level independently 
of B{}i^ —>■ T~^ u-r) , in agreement with the expected exclusion sensitivity of 57 GeV/c^. 



5 Conclusions 

The search for pair-produced charged Higgs bosons in the three final states T^UrT^Ur, 
csT~Ur and esse has been updated using 56.9 pb^^ of data collected at y/s = 181 — 184 GeV. 
No evidence of Higgs boson production was found and new mass limits were set as a function 
of B(H"'"— ^-r+i/,-). When combined with data recorded at centre-of-mass energies from 130 
to 172 GeV, charged Higgs bosons with masses below 59 GeV/c^ are excluded at 95% C.L. 
independently of B(H+^r+z/T-). 
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Figure 4: Limit at 95% C.L. on the mass of charged Higgs bosons as a function of B(H"'" ^t+z^t-) 
from each of the three individual analyses. Shown are the expected (dash-dotted curves) and 
observed (shaded areas) 95% C.L. exclusion domains. 
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Figure 5: Limit at 95% C.L. on the mass of charged Higgs bosons as a function of B(H'^ ^t+z^t-) 
from the combination of the three individual analyses. Shown are the expected (dash-dotted curve) 
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